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Abstract 

Multijet production rates in neutral current deep inelastic scattering have been 
measured in the range of exchanged boson virtualities 10 < Q 2 < 5000 GeV 2 . 
The data were taken at the ep collider HERA with centre-of-mass energy y/s = 
318 GeV using the ZEUS detector and correspond to an integrated luminosity of 
82.2 pb" 1 . Jets were identified in the Breit frame using the kr cluster algorithm 
in the longitudinally invariant inclusive mode. Measurements of differential dijet 
and trijet cross sections are presented as functions of jet transverse energy (-E 1 ^), 
pseudorapidity (?7lab) an d Q 2 with E^ B > 5 GeV and —1 < < 2.5. Next- 
to-leading-order QCD calculations describe the data well. The value of the strong 
coupling constant a s (Mz), determined from the ratio of the trijet to dijet cross 
sections, is a s (M z ) = 0.1179 ± 0.0013 (stat.)+°;°°^ (exp.)i°;°°^j (th.). 
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1 Introduction 



Measurements of multijet production from initial-state hadrons and leptons have been 
carried out previously in collisions at the SPS [1], the ISR [2], the TEVATRON [3] and 
at LEP [4] as well as in photoproduction [5] and deep inelastic scattering (DIS) [6] at 
HERA. Multijet production in DIS at HERA has been used to test the predictions of 
perturbative QCD (pQCD) calculations over a large range of four-momentum transfer 
squared, Q 2 [7]. Recently, the ZEUS and HI collaborations have determined the strong 
coupling constant a s from a variety of measurements of jet production and jet properties 
in both DIS [6,8-14] and photoproduction [15]. 

At leading order (LO) in a s , dijet production in neutral current DIS proceeds via the 
boson-gluon-fusion (BGF, V*g -> qq with V = 7, Z°) and QCD-Compton (QCDC, 
V*q — > qg) processes. Events with three jets can be seen as dijet processes with an 
additional gluon radiation or splitting of a gluon into a quark-antiquark pair and are 
directly sensitive to 0(a 2 ) QCD effects. The higher sensitivity to a s and the large number 
of degrees of freedom of the trijet final state allow detailed testing of QCD predictions. 

In the present analysis, the differential cross sections for the trijet production have been 
measured with high statistical precision. Measurements of the inclusive trijet cross section 
as a function of Q 2 and the jet transverse energy, E^ B , in the Breit frame and the jet 
pseudorapidity, ?7lab> m the laboratory frame are presented. Predictions of pQCD at 
next-to-leading order (NLO) are compared to the measurements. In addition, the analysis 
includes the first a s determination using the cross-section ratio of trijet to dijet production, 
R3/2, at HERA. In this ratio, correlated experimental and theoretical uncertainties cancel, 
allowing for an extension of the measurement to low Q 2 . 



2 Experimental set-up 

The data used in this analysis were collected during the 1998-2000 running period, when 
HERA operated with protons of energy E p = 920 GeV and electrons or positrons 1 of 
energy E e = 27.5 GeV, and correspond to an integrated luminosity of 82.2 ±1.9 pb" 1 . A 
detailed description of the ZEUS detector can be found elsewhere [16, 17]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are measured in the central tracking detector (CTD) [18], which oper- 
ates in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organised in nine superlayers covering the 

1 In the following, the term "electron" denotes generically both the electron (e~) and the positron (e + ). 
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polar-angle 2 region 15° < 9 < 164°. The transverse momentum resolution for full-length 
tracks can be parameterised as <j(pt)/pt = 0.0058pr © 0.0065 © 0.0014/pr, with pt in 
GeV. The tracking system was used to measure the interaction vertex with a typical 
resolution along (transverse to) the beam direction of 0.4 (0.1) cm and also to cross-check 
the energy scale of the calorimeter. 

The high-resolution uranium-scintillator calorimeter (CAL) [19] covers 99.7% of the total 
solid angle and consists of three parts: the forward (FCAL), the barrel (BCAL) and 
the rear (RCAL) calorimeters. Each part is subdivided transversely into towers and 
longitudinally into one electromagnetic section (EMC) and either one (in RCAL) or two (in 
BCAL and FCAL) hadronic sections (HAC). The smallest subdivision of the calorimeter 
is called a cell. Under test-beam conditions, the CAL single-particle relative energy 
resolutions were a(E)/E = 0.18 /VE for electrons and a(E)/E = 0.35/ \fE for hadrons, 
with E in GeV. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > e^p. The 
resulting small angle energetic photons were measured by the luminosity monitor [20], a 
lead-scintillator calorimeter placed in the HERA tunnel at Z = —107 m. 

3 Kinematics and event selection 

A three-level trigger system was used to select events online [17,21-23]. Neutral current 
DIS events were selected by requiring that the scattered electron with energy more than 

4 GeV was measured in the CAL [24]. 

The offline kinematic variables Q 2 (four- momentum transfer squared), x-q] (Bjorken scaling 
variable) and y = Q 2 / (sxbj) (s is the centre-of-mass energy squared) were reconstructed by 
the electron (e), double angle (DA) [25] and Jacquet-Blondel (JB) [26] methods. The angle 
of the hadronic system, 7h a d, corresponds, in the quark-parton model, to the direction of 
the scattered quark and was reconstructed from the CAL measurements of the hadronic 
final state. 

If 7had w as less than 90° and the scattered-electron track could be well reconstructed by 
the CTD, the DA method was used; otherwise, the electron method was used. The offline 
selection of DIS events was similar to that used in a previous ZEUS measurement [27] 
and was based on the following requirements: 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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• E' e >10 GeV, where E' e is the scattered-electron energy after correction for energy loss 
in inactive material in front of the CAL, to achieve a high-purity sample of DIS events; 

• y e <0.6, where y e is y reconstructed by the electron method, to reduce the photopro- 
duction background; 

• Z/jb>0.04, where j/jb is y reconstructed by the JB method, to ensure sufficient accuracy 
for the DA reconstruction of Q 2 ; 

• cos7had<0.7, to ensure good reconstruction of jets in the Breit frame; 

• 40<^ i (£' — Pz)i<Q0GeV, where the sum runs over all CAL energy deposits. The 
lower cut removed background from photoproduction and events with large initial- 
state QED radiation. The higher cut removed cosmic-ray background; 

• |Z vcrtcx |<50 cm, where Z vertex is the reconstructed primary vertex Z-position, to select 
events consistent with ep collisions; 

• |X|>13 or |V|>7 cm, where X and Y are the impact positions of the scattered electron 
on the RCAL, to avoid the low-acceptance region adjacent to the rear beampipe. 

The kinematic range of the analysis is defined as: 

10 < Q 2 < 5000 GeV 2 and 0.04 < y < 0.6. 

Jets were reconstructed using the hr cluster algorithm [28] in the longitudinally invariant 
inclusive mode [29]. The jet search was conducted in the Breit frame [30]. For each 
event, the jet search was performed using a combination of track and CAL information, 
excluding the cells and the track associated with the scattered electron. The selected 
tracks and CAL clusters were treated as massless Energy Flow Objects (EFOs) [31]. The 
clustering of objects was done according to the Snowmass convention [32]. 

The jet phase space is defined by selection cuts on the jet pseudorapidity 77lab m the 
laboratory frame and on the jet transverse energy E i ^ t B in the Breit frame: 

- 1 <^lab< 2 - 5 and E^ B >5GeV. 

Events with two (three) or more jets were selected by requiring the invariant mass of the 
two (three) highest E^ B jets to be: 

M 2jets(3jcts) > 25 GeV. 

These requirements were necessary to ensure a reliable prediction of the cross sections at 
NLO (see Section EJ). 

After all cuts, 37089 events with two or more jets (dijets) and 13665 events with three or 
more jets (trijets) remained. 
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4 Monte Carlo simulation 



Monte Carlo (MC) simulations were used to correct the data for detector effects, ineffi- 
ciencies of the event selection and that of the jet reconstruction, as well as for QED effects. 
Neutral current DIS events were generated using the Ariadne 4.08 program [33] and the 
Lepto 6.5 program [34] interfaced to Heracles 4.5.2 [35] via Django 6.2.4 [36]. The 
Heracles program includes QED effects up to C(a| M ). In case of Ariadne, the QCD 
cascade is simulated using the colour-dipole model [37], whereas for Lepto, the matrix 
elements plus parton shower model is used. Both models use the Lund string model [38], 
as implemented in Jetset 7.4 [39,40], for hadronisation. 

The ZEUS detector response was simulated with a program based on Geant 3.13 [41]. 
The generated events were passed through the detector simulation, subjected to the same 
trigger requirements as the data, and processed by the same reconstruction and offline 
programs. 

Measured distributions of kinematic variables are well described by both the Ariadne 
and Lepto MC models after reweighting in Q 2 . The Lepto simulation gives a better 
overall description of the E^ B and invariant mass distributions. Therefore, the events 
generated with the Lepto program were used to determine the acceptance corrections. 
The events generated with Ariadne were used to estimate the systematic uncertainty 
associated with the treatment of the parton shower. 

5 NLO QCD calculations 

The NLO calculations were carried out in the MS scheme for five massless quark fla- 
vors with the program Nlojet [42] using CTEQ6 [43], CTEQ4 [44], MRST99 [45] and 
ZEUS-S [46] for the proton parton density functions (PDFs). Nlojet allows a compu- 
tation of the trijet production cross sections to next-to-leading order, i.e. including all 
terms up to 0(a z s ). It was checked that the LO and NLO calculations from Nlojet agree 
with those of Disent [47] at the 1-2% level for the dijet cross sections [22,23]. 

For comparison with the data, the CTEQ6 parameterisations of the proton PDFs were 
used and the renormalisation and factorisation scales were both chosen to be (E^, + Q 2 )/4, 
where for dijets (trijets) Et is the average Et of the two (three) highest Et jets in a given 
event. The strong coupling constant was set to the value used in the CTEQ6 analysis, 
a s (Mz) = 0.1179, and evolved according to the two-loop solution of the renormalisation 
group equation. 

The NLO QCD predictions were corrected for hadronisation effects using a bin-by-bin 
procedure. Hadronisation correction factors were defined for each bin as the ratio of the 
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parton- to hadron-level cross sections and were calculated using the Lepto MC program. 
The correction factors Chad were typically in the range 1.15 — 1.35 for most of the phase 
space. 

6 Corrections and systematic uncertainties 

The jet transverse energy was corrected for energy losses in the inactive material in front 
of the CAL using the samples of MC simulated events [14]. The cross sections for jets of 
hadrons in bins of Q 2 , E^ B and ^lab were obtained by applying a bin-by-bin correction 
to the measured jet distributions using the Lepto program. The corrections take into 
account the efficiency of the trigger, the selection criteria and the purity and efficiency 
of the jet reconstruction. Additional corrections for QED effects, Cqed ; calculated using 
Heracles, were applied to the measured cross sections, aBom = o"moas. ■ Cqed- 

A detailed study of the sources contributing to the systematic uncertainties of the mea- 
surements was performed [22,23]. The main sources contributing to the systematic uncer- 
tainties are listed below (typical values of the systematic uncertainties in the dijet cross 
section and cross-section ratio R3/2 are indicated in parentheses): 

• jet-pseudorapidity cut- a change of ±0.1 (corresponding to the resolution) in the ?7lab 
cuts imposed on the jets in the laboratory frame for both data and MC simulated 
events (1%,1%); 

• jet transverse energy and invariant mass cuts - E^ B and M2j e t s (Af3j ets ) were simulta- 
neously varied by the corresponding resolution near the cuts for both data and MC 
simulated events. Along with the previous systematic check, this takes into account the 
effect of the remaining differences between the data and the MC simulation (3%, 3%); 

• use of different parton shower model - using Ariadne instead of Lepto to evaluate 
the acceptance corrections (2%, 4%); 

• the absolute energy scale of the CAL- varying E^ B by its uncertainty of ±1%(>10 GeV) 
and ±3%(<10GeV) for MC events [10] (6%,3.5%). 

The systematic uncertainties not associated with the absolute energy scale of the CAL 
were added in quadrature to the statistical uncertainties and are shown on the figures 
as error bars. The uncertainty due to the absolute energy scale of the CAL is highly 
correlated from bin-to-bin and is shown separately as a shaded band. The total system- 
atic uncertainty and the uncertainty due to the absolute energy scale are also shown in 
Tables [U-ini 

The main contributions to the theoretical uncertainties of the NLO QCD predictions 
are: 
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• uncertainties in the proton PDFs, which were estimated by repeating the calculations 
using 40 additional sets obtained under different theoretical assumptions as part of 
the CTEQ6 release (2.5%,2%); 

• uncertainties in the correction factors, Chad, which were estimated by using the Ari- 
adne program instead of LEPTO (6%, 4%); 

• uncertainties due to terms beyond NLO, which were estimated by varying both ixr and 
Hf between (E 2 + Q 2 ) and (E 2 + Q 2 )/16 (10%,7%). 

The total theoretical uncertainty was obtained by adding in quadrature the individual 
uncertainties listed above. 

7 Results 

7.1 Differential cross sections 

The differential trijet cross sections as functions of E^ B are presented in Fig. ^ and in 
Tables The three highest E^ B jets were ordered in E^ B (E^ t B 1 >E^ t B >E^ t B ). The 
observed decrease of the cross section for the first jet towards small values of Ej! B is 
caused by the Ej, B ordering combined with the requirement that the second and third 
jet have E^ B >5 GeV. For the second jet, a similar but less pronounced effect is observed. 
The NLO predictions using Nlojet, corrected for hadronisation effects, are compared to 
the data in Fig. The QCD predictions provide a good description of both the shape 
and magnitude of the measured cross sections, even at low E^ B . 

Figure El and Tables EHHl show the differential trijet cross sections as functions of ^lab- 
The three highest E^ B jets were ordered in j^ab (^lab^lab^lab)- Figure El and 
Tables Eland [HI show both the differential dijet and trijet cross section as functions of Q 2 . 
In Figs. |U and El the data are generally well described by the NLO QCD predictions. The 
largest difference is a slightly different slope of the t/lab dependence of the third jet. 

7.2 Cross-section ratio and determination of a s 

Figure 13] and Table |H1 show the cross-section ratio R3/2 of the trijet cross section to the 
dijet cross section, as a function of Q 2 . The correlated systematic and the renormalisation 
scale uncertainties largely cancel in the ratio. The agreement between the data and 
NLO predictions is good. The total experimental and theoretical uncertainties are about 
5% and 7%, respectively. These uncertainties are substantially reduced with respect to 
those of the di- and trijet cross sections. In particular, at low Q 2 (Q 2 < 100 GeV 2 ), the 
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theoretical uncertainties are reduced by as much as a factor of four. This reduction allows 
the determination of a s (Mz) at a much lower Q 2 than in previous analyses [10, 14]. 

The measurement of R3/2 as a function of Q 2 was used to determine a s (Mz) with a 
method similar to that of a previous ZEUS publication [9]: 

• the NLO QCD calculation of R3/2 was performed using the five sets of proton PDFs of 
the CTEQ4 A-series [44] 3 . The value of a s (Mz) used in each partonic cross-section 
calculation was that associated with the corresponding set of PDFs: 0.110, 0.113, 
0.116, 0.119, 0.122; 

• for each bin, i, in Q 2 , the NLO QCD calculations, corrected for hadronisation effects, 
were used to parameterise the a s (Mz) dependence of R3/2 according to the functional 
form: 

[R 3/ 2(a s (M z ))} 1 = C\ ■ a s {M z ) + C\ ■ a 2 (M z ), (1) 

where C\ and C l 2 are fitting parameters. This simple parameterisation gives a good 
description of the a s (M z ) dependence of R^^iQ 2 ) over the entire a s range spanned 
by the PDF sets; 

• a value of a s (Mz) was then determined in each bin of Q 2 , as well as in the entire 
Q 2 region, by a x 2 -fit of the measured R^^iQ 2 ) values using the parameterisation in 
Eq. (0). 

This procedure correctly handles the complete a s -dependence of the NLO differential cross 
sections (the explicit dependence coming from the partonic cross sections and the implicit 
dependence coming from the PDFs) in the fit, while preserving the correlation between 
a s and the PDFs. Taking into account only the statistical uncertainties on the measured 
cross-section ratio, a s (M z ) is determined to be a s (M z ) = 0.1179±0.0013(stat.). 

Figure shows the sensitivity of the cross-section ratio R 3 / 2 to the value of a s . Figure 
and Table ITUl show a s [Mz) determined in the five bins of Q 2 . 

As a cross-check of the extracted value of a s (M z ), the fit procedure was repeated by using 
the three sets of the MRST99 PDF corresponding to a s {Mz) equal to 0.1125, 0.1175 and 
0.1225. The result is a s {M z ) = 0.1178±0.0010(stat.) 

In addition, the NLO QCD analysis used to obtain the ZEUS-S PDF [46] was re- 
peated to obtain a set of five PDFs corresponding to the values of a s (Mz). 0.115, 0.117, 
0.119, 0.121, 0.123. These sets were used in the current analysis yielding a s (M z ) = 
0.1191±0.0010(stat.), in good agreement with the other determinations. 

3 The CTEQ4 PDF was chosen because the CTEQ6 does not provide PDF sets obtained with different 
a s (Mz) values and therefore cannot be used for the determination of a s - 
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The experimental and theoretical uncertainties of the extracted value of a s (M z ) were 
evaluated by repeating the analysis above for each systematic check, as described in 
Section IH1 The main contributions to the experimental systematic uncertainty were: 

• jet pseudorapidity cut (1} 5%); 

• jet transverse energy and invariant mass cuts (^2% % ); 

• use of different parton shower model (—2%) ; 

• the absolute energy scale of the CAL (^2^%)- 

The main contributions to the theoretical uncertainty are: 

• uncertainties in the proton PDFs {^} 2 %°)'i 

• uncertainties in the correction factor, Chad (+2%); 

• uncertainties due to terms beyond NLO 

The value of a s (M z ) as determined from the measurements of R3/2 is therefore: 
a s {M z ) = 0.1179 ±0.0013 (stat.)t° S (e«p.)±aoo8 (th.). 

The result is in good agreement with recent determinations at HERA [6,8-15] and the 
current world average of a s (M z ) = 0.1182±0.0027 [48]. 



8 Summary 

Differential dijet and trijet cross sections have been measured with high precision in 
neutral current deep inelastic scattering for 10 < Q 2 < 5000 GeV 2 at HERA using the 
ZEUS detector. The inclusive trijet cross section has been measured as a function of E^ B , 
77lab an d Q 2 ■ The ratio R3/2 of the trijet and dijet cross sections has been measured as 
a function of Q 2 . The predictions of perturbative QCD calculations in next-to-leading 
order give a good description of the dijet and trijet cross sections and the cross-section 
ratio R3/2 over the whole range of Q 2 . The cancellation of uncertainties in the ratio, 
in particular those from theory, allow the extraction of a s with good precision down 
to Q 2 of 10 GeV 2 . The value of the strong coupling constant a s was measured to be 
a s [M z ) = 0.1179 ± 0.0013 (stat.)lHol (exp-)±aoo46 ( th -)> in g° od agreement with the 
current world average value and previous determinations of a s (M z ) at HERA. 



9 Acknowlegements 

We thank the DESY Directorate for their strong support and encouragement. The re- 
markable achievements of the HERA machine group were essential for the successful 



8 



completion of this work and are greatly appreciated. The design, construction and instal- 
lation of the ZEUS detector has been made possible by the effort of many people who are 
not listed as authors. We would like to thank Z. Nagy for useful discussions. 



9 



References 

[1] UAl Coll., G. Arnison et al., Phys. Lett. B 158, 494 (1985); 
UA2 Coll., J. A. Appel et al, Z. Phys. C 30, 341 (1986). 

[2] CMOR Coll., A. L. S. Angelis, Nucl. Phys. B 303, 569 (1988); 

Axial Field Spectrometer Coll., T. Akesson et al, Z. Phys. C 32, 317 (1986). 

[3] DO Coll., B. Abbott et al., Phys. Rev. Lett. 86, 1955 (2001); 
DO Coll., S. Abachi et al, Phys. Rev. D 53, 6000 (1996). 

[4] P. Lutz, Proc. Les Arcs 1996, QCD and high energy hadronic interactions, 
pp. 223-228. (1996). Prepared for 31st Rencontres de Moriond: QCD and 
High-Energy Hadronic Interactions, Les Arcs, France, 23-30 Mar 1996; 
Z. Nagy and Z. Trocsanyi, Nucl. Phys. Proc. Suppl. 74, 44 (1999). 

[5] ZEUS Coll., J. Breitweg et al, Phys. Lett. B 443, 394 (1998). 

[6] HI Coll., C. Adloff et al., Eur. Phys. J. C 6, 575 (1999). 

[7] HI Coll., C. Adloff et al., Phys. Lett. B 515, 17 (2001). 

[8] ZEUS Coll., M. Derrick et al., Phys. Lett. B 363, 201 (1995). 

[9] ZEUS Coll., J. Breitweg et al., Phys. Lett. B 507, 70 (2001). 

[10] ZEUS Coll., S. Chekanov et al, Phys. Lett. B 547, 164 (2002). 

[11] ZEUS Coll., S. Chekanov et al., Phys. Lett. B 558, 41 (2003). 

[12] HI Coll., C. Adloff et al., Eur. Phys. J. C 19, 289 (2001). 

[13] HI Coll., C. Adloff et al., Eur. Phys. J. C 21, 33 (2001). 

[14] ZEUS Coll., S. Chekanov et al., Nucl. Phys. B 700, 3 (2004). 

[15] ZEUS Coll., S. Chekanov et al., Phys. Lett. B 560, 7 (2003). 

[16] ZEUS Coll., M. Derrick et al., Phys. Lett. B 293, 465 (1992). 

[17] ZEUS Coll., U. Holm (ed.), The ZEUS Detector. Status Report (unpublished), 

DESY (1993), available on http://www-zeus.desy.de/bluebook/bluebook.html. 

[18] N. Harnew et al., Nucl. Inst. Meth. A 279, 290 (1989); 

B. Foster et al, Nucl. Phys. Proc. Suppl. B 32, 181 (1993); 
B. Foster et al, Nucl. Inst. Meth. A 338, 254 (1994). 

[19] M. Derrick et al, Nucl. Inst. Meth. A 309, 77 (1991); 
A. Andresen et al., Nucl. Inst. Meth. A 309, 101 (1991); 
A. Caldwell et al, Nucl. Inst. Meth. A 321, 356 (1992); 
A. Bernstein et al, Nucl. Inst. Meth. A 336, 23 (1993). 



10 



[20] J. Andruszkow et al, Preprint DESY-92-066, DESY, 1992; 
ZEUS Coll., M. Derrick et al, Z. Phys. C 63, 391 (1994); 
J. Andruszkow et al, Acta Phys. Pol. B 32, 2025 (2001). 

[21] ZEUS Coll., J. Breitweg et al., Z. Phys. C 74, 207 (1997). 

[22] N. Krumnack, Ph.D. Thesis, University of Hamburg, 2004. 

[23] L. Li, Ph.D. Thesis, University of Wisconsin-Madison, 2005. 

[24] ZEUS Coll., J. Breitweg et al., Eur. Phys. J. C 11, 427 (1999). 

[25] S. Bentvelsen, J. Engelen and P. Kooijman, Proc. Workshop on Physics at HERA, 
W. Buchmiiller and G. Ingelman (eds.), Vol. 1, p. 23. Hamburg, Germany, DESY 
(1992). 

[26] F. Jacquet and A. Blondel, Proceedings of the Study for an ep Facility for Europe, 
U. Amaldi (ed.), p. 391. Hamburg, Germany (1979). Also in preprint DESY 79/48. 

[27] ZEUS Coll., S. Chekanov et al, Eur. Phys. J. C 23, 13 (2002). 

[28] S.Catani et al., Nucl. Phys. B406, 187 (1993). 

[29] S.D. Ellis and D.E. Soper, Phys. Rev. D 48, 3160 (1993). 

[30] R.P. Feynman, Photon-Hadron Interactions. Benjamin, New York, 1972; 
K.H. Streng, T.F. Walsh, P.M. Zerwas, Z. Phys. C 2, 237 (1979). 

[31] ZEUS Coll., S. Chekanov et al, Preprint hep-ex/0404033, 2004. 

[32] J.E. Huth et al., Research Directions for the Decade. Proceedings of Summer Study 
on High Energy Physics, 1990, E.L. Berger (ed.), p. 134. World Scientific (1992). 
Also in preprint FERMILAB-CONF-90-249-E. 

[33] L. Ldnnblad, Comp. Phys. Comm. 71, 15 (1992). 

[34] G. Ingelman, A. Edin and J. Rathsman, Comp. Phys. Comm. 101, 108 (1997). 

[35] A. Kwiatkowski, H. Spiesberger and H.-J. Mohring, Comp. Phys. Comm. 

69, 155 (1992). Also in Proc. Workshop Physics at HERA, 1991, DESY, Hamburg. 

[36] K. Charchula, G.A. Schuler and H. Spiesberger, Comp. Phys. Comm. 
81, 381 (1994). 

[37] G. Gustafson and U. Pettersson, Nucl. Phys. B 306, 746 (1988). 

[38] B. Andersson et al., Phys. Rep. 97, 31 (1983). 

[39] M. Bengtsson and T. Sjostrand, Comp. Phys. Comm. 46, 43 (1987). 

[40] T. Sjostrand, Comp. Phys. Comm. 82, 74 (1994). 

[41] R. Brun et al, geant3, Technical Report CERN-DD/EE/84-1, CERN, 1987. 



11 



[42] Z. Nagy and Z. Trocsanyi, Phys. Rev. Lett. 87, 082001 (2001). 

[43] J. Pumplin et al., Preprint hep-ph/0201195, 2002. 

[44] H.L. Lai et al, Phys. Rev. D 55, 1280 (1997). 

[45] A.D. Martin et al, Eur. Phys. J. C 14, 133 (2000). 

[46] ZEUS Coll., S. Chekanov et al., Phys. Rev. D 67, 012007 (2003). 

[47] S. Catani and M. H. Seymour, Nucl. Phys. B 485, 291 (1997). 

[48] S. Bethke, J. Phys. G 26, R27 (2000); 

Updated in S. Bethke, Preprint hep-ex/0407021, 2004. Talk given at 7th DESY 
Workshop on Elementary Particle Theory: Loops and Legs in Quantum Field 
Theory, Zinnowitz, Germany, 25-30 Apr 2004. 



12 



T,B 

(GeV) 


(pb/GeV) 


5 - 8 


5.22 


8 - 12 


20.2 


12 - 16 


14.4 


16 - 20 


6.05 


20 - 25 


2.19 


25 - 30 


0.828 


30 - 40 


0.222 


40 - 60 


0.047 



± 0.29 
± 0.5 
± 0.4 
± 0.22 
± 0.11 
± 0.068 
± 0.027 
± 0.012 



Jsyst 



+0.75 
-0.89 
+3.0 
-3.5 
+1.6 
-1.4 
+0.68 
-0.54 
+0.32 
-0.21 
+0.12 
-0.10 
+0.036 
-0.026 
+0.004 
-0.012 



+0.63 
-0.50 
+3.0 
-2.9 
+1.6 
-1.3 
+0.42 
-0.44 

+0.20 
-0.16 
+0.063 
-0.065 
+0.021 
-0.017 
+0.001 
-0.004 



Cqed 

0.97 
0.96 
0.96 
0.96 
0.96 
0.97 
0.97 
1.07 



Chad 
2.11 

1.33 
1.21 
1.19 
1.19 
1.22 
1.20 
1.31 



Table 1: The inclusive trijet cross section da/dE^g for jets of hadrons in the 
Breit frame, selected using the k? cluster algorithm in the longitudinally invariant 
inclusive mode. The statistical (S sta t), systematic (S syst ) and the absolute energy 
scale uncertainties (Ses) ar ^ shown separately. The multiplicative correction factors 
for QED radiative effects (Cqed), applied to the data, and for hadronisation effects 
(Chadj? applied to the NLO predictions, are shown in the last two columns. 



77ijct,2 
T.B 

(GeV) 


da/dE^g <5 sta t 4 ys t <5es 
(pb/GeV) 


Cqed 


Chad 


5 - 8 
8 - 12 
12 - 16 
16 - 20 
20 - 25 
25 - 30 
30 - 40 


24.0 ± 0.6 tH ±H 
20.7 ±0.4 j£° t|° 
6.70 ±0.24 jftjjjj 
oifi +013 + - 14 + - 13 

Z.IO IC U.IO _ 2 2 -0.14 

0.780 ± 0.067 ^;?o 7 -oZ 
225 ± 038 +0 ' 034 +a005 

U.ZZd IE U.UOO „ .005 -0.005 

n 0^)4 ±0014 + - 013 +a004 
u.ud^t in u.ui^t _ .006 -o.ooi 


0.96 
0.97 
0.96 
0.95 
0.97 
0.97 
1.02 


1.48 
1.22 
1.28 
1.32 
1.34 
1.40 
1.51 



Table 2: The inclusive trijet cross section da/dE^g for jets of hadrons in the 
Breit frame, selected using the kr cluster algorithm in the longitudinally invariant 
inclusive mode. Other details are as in the caption to TableU\ 
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T,B 

(GeV) 


aa / dhi TB d st3u t o sys t oes 
(pb/GeV) 


Cqed 


Chad 


5 - 8 
8 - 12 
12 - 16 
16 - 20 
20 - 25 


52.8 ±0.8 tH 
8.06 ± 0.25 t\i t\\ 

i.ii ±0.09 t ^ ±85S 

0.208 ± 0.039 

0.052 ± 0.020 ±g:gs t£8E 


0.96 
0.96 
0.97 
0.98 
0.94 


1.28 

1.49 
1.66 
1.77 
1.77 



Table 3: The inclusive trijet cross section da/dE^g for jets of hadrons in the 
Breit frame, selected using the kx cluster algorithm in the longitudinally invariant 
inclusive mode. Other details are as in the caption to Tabled 



jet.l 
^LAB 


da/drfl\ B 5 stat S syst S E s 
(Pb) 


Cqed 


Chad 


0.5 - 1.0 
1.0 - 1.5 
1.5 - 2.0 
2.0 - 2.5 


19.1 ±1.2 tH 

61.2 ± 2.0 tl 3 7 tH 

i ok _i_ o +15 +15 
1ZO 3= O _ 2 q _ 13 

186 ±4 tH tf 8 


0.95 
0.96 
0.96 
0.96 


2.18 

1.51 
1.31 
1.17 



Table 4: The inclusive trijet cross section da/drf^ B for jets of hadrons in the 
Breit frame, selected using the kx cluster algorithm in the longitudinally invariant 
inclusive mode. Other details are as in the caption to Tabled 



jet,2 
?7lAB 


dcr/drf^l 
(pb) 


Ostat 


^syst 


5es 


Cqed 


Chad 


-1.0 - -0.5 


0.34 


± 0.14 


+0.04 
-0.18 


+0.03 
-0.09 


0.91 


8.72 


-0.5 - 0.0 


13.6 


± 1.0 


+3.0 
-2.4 


+2.9 
-2.1 


0.98 


2.30 


0.0 - 0.5 


62.1 


± 2.1 


+9.3 
-11.0 


+9.1 
-8.4 


0.96 


1.50 


0.5 - 1.0 


115 


± 3 


+14 
-14 


+14 
-13 


0.97 


1.29 


1.0 - 1.5 


108 


± 3 


+12 
-14 


+ 11 
-10 


0.96 


1.22 


1.5 - 2.0 


73.3 


± 2.2 


+7.3 
-7.8 


+7.3 
-6.3 


0.96 


1.19 


2.0 - 2.5 


20.3 


± 1.4 


+5.6 
-2.1 


+1.7 
-1.6 


0.97 


1.16 



Table 5: The inclusive trijet cross section da/drf^ AB for jets of hadrons in the 
Breit frame, selected using the kx cluster algorithm in the longitudinally invariant 
inclusive mode. Other details are as in the caption to Tabled 
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jet, 3 
^LAB 


da/drflU 
(Pb) 


"stat 


-1.0 - -0.5 


63.7 


± 2.3 


-0.5 - 0.0 


112 


± 3 


0.0 - 0.5 


108 


± 3 


0.5 - 1.0 


68.6 


± 2.1 


1.0 - 1.5 


31.4 


± 1.5 


1.5 - 2.0 


6.21 


± 0.67 



-'syst 



+9.2 
-11.0 

+15 
-14 

+12 
-15 
+8.7 
-7.6 
+2.3 
-4.6 
+1.50 
-0.67 



+9.0 
-9.3 

+15 
-12 

+12 
-11 
+7.4 
-6.5 
+2.1 
-2.7 
+0.61 
-0.61 



Cqed 

0.96 
0.96 
0.96 
0.97 
0.94 
0.93 



Chad 
1.32 

1.24 
1.25 
1.32 
1.44 
1.50 



Table 6: The inclusive trijet cross section da/drf^^ for jets of hadrons in the 
Breit frame, selected using the kr cluster algorithm in the longitudinally invariant 
inclusive mode. Other details are as in the caption to Tabled 



Q 2 

(GeV 2 ) 


(da/dQ 2 )dijet ^stat ^syst °~ES 

(pb/GeV 2 ) 


Cqed 


Chad 


10 - 35 
35 - 85 
85 - 220 
220 - 700 
700 - 5000 


9.70 ±0.10 ±lf 4 t° f e 
2.50 ±0.04 US t° f 8 
0.649 ±0.011 ±g:gJS iaoli 
nin4 + 0002 +0 - 006 +a006 

V.lVi 3= U.UUZ -0.006 -0.006 

nnn4 no +Q 000 1 4 + - 00015 +0-00014 
U.UU^UO 33 U.UUUl^ -0.00024 -0.00014 


0.97 
0.95 
0.96 
0.94 
0.92 


1.18 
1.16 
1.15 
1.12 
1.09 



Table 7: The inclusive dijet cross section (da / 'dQ 2 )dij C t for jets of hadrons in the 
Breit frame, selected using the k? cluster algorithm in the longitudinally invariant 
inclusive mode. Other details are as in the caption to Tabled 



Q 2 

(GeV 2 ) 


(da /dQ 2 ) trijet 5 st at #syst <^ES 

(pb/GeV 2 ) 


Cqed 


Chad 


10 - 35 
35 - 85 
85 - 220 
220 - 700 
700 - 5000 


3.94 ±0.08 t° H t° H 
0.94 ±0.02 I} 2 lg£ 
0.227 ± 0.007 iffi %™ 
0320 ± 001 3 +0 0036 + 0027 

U.UOZU 33 U.UUIO -0.0031 -0.0026 

00112 ± 00007 +00000 8 +0.00008 

U.UU11Z 33 U.UUUU/ -0.00014 -0.00009 


0.98 
0.95 
0.96 
0.94 
0.92 


1.35 

1.31 
1.32 
1.35 
1.33 



Table 8: The inclusive trijet cross section (da/dQ 2 ) trijet for jets of hadrons in the 
Breit frame, selected using the k? cluster algorithm in the longitudinally invariant 
inclusive mode. Other details are as in the caption to TableU\ 
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Q 2 

(GeV 2 ) 


R3/2 "stat "syst "ES 


Cqed 


Chad 


10 - 35 
35 - 85 
85 - 220 
220 - 700 
700 - 5000 


406 ± 008 +0 ' 014 + - 013 

U.^UU 33 U.UUO _ .027 -0.017 

n o 7 p- _un 010 +0 ' 014 +0 - 014 

U.O/O 33 U.UIU _ .019 -0.010 
U.OOU 33 U.U11 _ .031 -0.015 

0.306 ±0.013 tffi 

0.279 ±0.018 + ^29 -ao?4 


1.01 
0.99 
1.00 
1.01 
1.01 


1.14 
1.02 
1.15 
1.21 
1.21 



Table 9: The ratio of inclusive trijet to dijet cross sections for jets of hadrons 
in the Breit frame, selected using the kx cluster algorithm in the longitudinally 
invariant inclusive mode. Other details are as in the caption to TableU\ 



Q 2 

(GeV 2 ) 


a s (M z ) 5 s tat 4x P <Wo 


10 - 35 
35 - 85 
85 - 220 
220 - 700 
700 - 5000 


0.1210 ± 0.0022 + -°oZt 
0.1148 ± 0.0024 lg£» tffi 
n 1 1 78 -4- n nno7 +0.0027 +0.0064 

U.llrtS ± 0.002^ _ .oo63 -0.0016 

n 1 171 + n nn39 +0-0043 +0.0068 

U.1I/1 zSz U.UUOJ _ .0056 -0.0023 

1 1 70 ± 0064 + - 0045 +°- 011 

U.ll/U 33 U.UUU4: -0.0064 -0.0037 


10 - 5000 


mi 7Q 4-nnm^ +0.0028 +0.0064 

U.llCy ± U.UUlJ -Q.0046 -0.0046 



Table 10: The a s (Mz) values as determined in this analysis. The statistical 
(S s ta.t), experimental (S exp ) and theoretical (S t heo) systematic uncertainties are shown 
separately. 
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Figure 1: (a) The inclusive trijet cross sections as functions of E^ B with the 

jets ordered in Ej, B . The cross sections of the second and third jet were scaled 
by the factors shown for readability. The inner error bars represent the statisti- 
cal uncertainties. The outer error bars represent the quadratic sum of statistical 
and systematic uncertainties not associated with the calorimeter energy scale. The 
shaded band indicates the calorimeter energy scale uncertainty. The predictions of 
perturbative QCD in next-to-leading order, corrected for hadronisation effects and 
using the CTEQ6 parameterisations of the proton PDFs, are compared to the data, 
(b), (c) and (d) show the ratio of the data to the predictions. The hatched band 
represents the renormalisation scale uncertainty of the QCD calculation. 
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Figure 2: (a) The inclusive trijet cross sections as functions ofrf^ AB with the jets 
ordered in //lab ■ ^ e cross sections of the second and third jet were scaled up for 
readability only. The predictions of perturbative QCD in next-to-leading order are 
compared to the data, (b), (c) and (d) show the ratio of the data to the predictions. 
Other details are as in the caption to Fig. QJ 
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Figure 3: (a) The inclusive dijet and trijet cross sections as functions of Q 2 . 
The predictions of perturbative QCD in next-to-leading order are compared to the 
data, (b) and (c) show the ratio of the data to the predictions. Other details are 
as in the caption to Fig. 
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Figure 4: (a) The ratio of inclusive trijet to dijet cross sections as a function of 
Q 2 . The predictions of perturbative QCD in next-to-leading order are compared to 
the data, (b) shows the ratio of the data to the predictions. Other details are as in 
the caption to Fig. 
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Figure 5: (a) The ratio of inclusive trijet to dijet cross sections as a function of 
Q 2 . The predictions of perturbative QCD in next-to-leading order using five sets of 
CTEQ4 PDF are compared to the data, (b) shows the a s [Mz) values determined 
from the ratio of inclusive trijet to dijet cross sections in different regions of Q 2 . 
The shaded band indicates the current world average value of a s {Mz) ■ The inner 
error bars represent the statistical uncertainty of the data. The outer error bars 
show the statistical and systematic uncertainties added in quadrature. The dashed 
error bars display the theoretical uncertainties. 
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